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Abstract 



00 ' The next generation of particle edm searches will be at such a high sensitivity that it will be 

o ■ 

^ ' possible for the results to be contaminated by a systematic error resulting from the interaction of 

^ ■ the motional {E x v/c) magnetic field with stray field gradients. In this paper we extend previous 

H : 

, work to present an analytic form for the frequency shift in the case of a rectangular storage vessel 

and discuss the implications of the result for the neutron edm experiment which will be installed 
at the SNS (Spallation Neutron Source) by the LANL collaboration 
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I. INTRODUCTION 

Searches for particle electric dipole moments (edm) are considered to be one of the most 
promising places to search for physics beyond the standard model. Current experiments 
have reached the sensitivity where they have to take into account a systematic effect due 
to the influence on the particle's magnetic dipole moment of a combination of the motional 
magnetic fleld, (^B^ = E xl^/c^ due to the motion of the particle in the applied static 
electric field, and gradients in the ambient magnetic field. For slow motion (adiabatic limit) 
the effect can be described as a geometric phase effect [1] whereas a description in terms 
of the Bloch-Siegert shift is valid for fast motions as well ([2]). A treatment vahd also for 
intermediate motions for particles moving in a cylindrical container has been given in [3], 
[4]. In this note we discuss some additional symmetries of the effect and present an analytic 
solution for the case of particles moving in a rectangular vessel. 
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II. SYMMETRY OF THE EFFECT 



It is easy to see [2], that the effect depends only on motion in a plane perpendicular to 
the direction of E (parallel motion produces no -B^). Equations (23 and 36) in [3] show- 
that in general the frequency shift linear in E which leads to the systematic error is given 
in terms of the spectrum of the velocity autocorrelation function by 



7^^ 



2c 
with 



Sy (uo) + 5'x (uJo) 



dy dx 



(1) 



/■oo 

Si {uJo) = / cos LUoTRi (t) dr (2) 
Jo 

Ri (r) = 2 / (x) dx (3) 
Jo 

ilJi{x)^{vi{t)vi{t-x)) (4) 



oo 



S, (ujo) = 2 / r^^dw (5) 

where cUg — 7-Bo is the Larmor frequency, 7 is the gyromagnetic ratio, Eg is the homogeneous 
magnetic field taken as parallel to E in the z direction and (cu) is defined by 

/oo 
cos Lox^i (uj) dw (6) 
-00 

. If the field is cylindrically symmetric 

dy dx 2 dz 

^^E-^l^iSyM + sAu^o)) (8) 

and the result depends only on for all frequencies and geometries of the orbits. In the 
case that the trajectories have cylindrical symmetry 

Sy (ujo) = S.^ (ujo) = S (ujo) (9) 

This symmetry will hold in the high frequency limit {cuo » cUr) which is determined by 
the short time behavior of the correlation function, and the result that the shift in this case 
depends only on ^ has been obtained in [2] , section IV B and by a different method in [4] 
section IV A. The advantage of the present treatment is obvious. 
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An interesting case that can arise in practice is the symmetry 



dz 



0, 

dx 



dx 



dBy 

dy 



» 



(11) 

(12) 



In the high frequency hmit this difference in the spectra will approach zero. 



III. MOTION IN A RECTANGULAR BOX 



A. Single velocity 

In a rectangular box with specular reflecting walls the orbits are straight hues reflecting 
at equal angles to the normal when they encounter a wall. With perpendicular walls parallel 
to the X and y axes the motion in each dimension will be independent of the other dimension 
and the angle with the normal will be preserved independently for the walls along x and 
y. At a wall collision the magnitude of velocity and the component parallel to the wall are 
unchanged and the perpendicular component of velocity changes sign. Thus the correlation 
function for a given velcocity component {i — x, y) starts at and after I collisions is 

^lJ^^r) = {v,{T)v,m=vU-l^ for IT, < r < {I + l)T, (13) 

with the time between collisions — L^/vi, being constant for each orbit and the length 
of the box in direction i. 

Thus the velocity correlation function for each direction is a square wave with the switch- 
ing points being equally spaced but whose exact timing depends on the distance of the 
starting point of a given orbit from the first wall collision. Averaging over these starting 
points proceeds as in equations 14-17 of [4]. In fact the orbits are exactly those characterized 
by q; = 7r/2 in that reference (see flg.2 [4]) and the results of that paper can be apphed to 
the present case by substituting 

a = 7i/2 (14) 
R = Li/2 (15) 



Thus equ. 42 [4] becomes 



4 



E 



- (tt (m + 1/2))^) 
-oo + 1/2))' |_ - (tt (m + 1/2))')' + 



(16) 
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with uj'l — uJoLi/2vi and = -^i/2A, A being the mean free path between gas coUisions. 
Introducing uj'^ — ujoLi/vi we have 



TT (m + 1/2))' 1^ (^(^/2 _ (2m + 1))')' + u'^^l^ 



(17) 



with lo — This is plotted in figure 1), which is to be compared with figure 3) of [4] for 
the case of a cyhnder. 

Equation (17) and figure 1) apply to motion in one dimension. For the full 2 dimensional 
problem it is necessary to add together suitably normalized forms of the function for each 
dimension. 

The value of Si {lOo = 0) = .083L^ agrees with that predicted by the diffusion theory for 

the single dimension contribution to a rectangular box (eqation 82, [3]), 

8 1 
-T } ^ = .083 



m=l,6,o 



(18) 



Equation 17 shows that Si {ujo) becomes independent of Lj as the frequency increases and 
the larger the damping (larger Zq), the higher the frequency where this occurs. 



B. Frequency shift averaged over Msixwell velocity distribution, the case of co- 
magnet omet er s . 

Both the neutron edm experiment carried out by Baker et al, [5] and that being developed 
by the Los Alamos collaboration [6], make use of co-magnctomctcrs, that is a gas of atoms 
occupying the same region as the ultra-cold neutrons and satisfying the Maxwell-Bolzmann 
velocity distribution. In this case we write the velocity as 

Vi = yf3{T) (19) 

with (3 (T) — -y/^ the most probable velocity in a volume. Then our one dimensional 
velocity Vi has the probability distribution 

P{y)dy 



The spectral function of the frequency shift can then be rewritten 



Si {uio, y, Q = Ll ^ 



1 



- (7r(m + l/2))^ 



{uf-{'K{2m+l)Yy^)y^ 
({uf-{7:{2m + l)fy^f + uf {llf) 



(20) 



(21) 
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where a;* = uJoLi/(5{T) , and Z* = Li/ {(5 {T)tc{T)) and we have speciahzed to the case 
where the gas coUision time, Tc is independent of velocity. This is vahd for the common case 
where the scattering cross section satisfies cr ~ l/f and holds in particular for the case of 
He^ diffusing in superfiuid He^ which is the co- magnetometer in the LANL experiment [6] , 
[7]. 



1. He^ colliding with phonons in superfiuid He^, the co-magnetometer in the LANL search 
for a neutron electric dipole moment. 

Since the velocity of the He^ is much less than the phonon velocity, the collision rate of 
the phonons with the He^ will be independent of velocity and the mean free path satisfies 
A = VTc- We obtain Tc from 

(T) = 3D (T) / {v% (22) 

with {v^)rp the mean square velocity in a volume of gas and D (T) has been measured at 
temperatures of interest [8] 

D (T) = ^cm' / sec (23) 

We now average the frequency shift over the Maxwell- Boltzman distribution for velocity in 
one dimension 

VI/, (a;:, T) = J e-y'S. {u^ Q dy (24) 

The results are plotted in figure 2) : 

The same result is shown as a function of temperature for fixed (normalized) frequency 
in fgure 3. (Note the frequency normalization is temperature dependent. 



C. A rectangular box with the two sides significantly different, 

Equations 17 and 24 refer to a single dimension. If the second dimension of the box has 
a length Ly, then equation 4 can be written {e — Ly/L^) , keeping the same normalization 
for Sy and lo*: 



Sy {ujo,y,ll) 
LI 



— £ 



oc 
m=—oo 



{u;f-{7r{2m + l)fy') y 



(^{e^uf - (tt (2m + 1))' y^^ + i^of 



(25) 
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Averaging this as in equation 24 we obtain, {"^y {uj*, T, I*) — {suj*, T, l^e)) 

* T) = T) + eHy {ul T) (26) 

Figure 4) shows the two dimensional result for e = 0.2, normalized to as a function of 
frequency normalized to for various temperatures while figure 5) shows the normalized 
shift vs temperature for various normalized frequencies. 

In figure 6) we show the contribution to the frequency shift for the two directions indepen- 
dently, compared to the results of numerical simulations, as well as the total frequency shift 
given by their sum for T=0.4K. The contribution of the short side (y) has been normalized 
to the long side (x) for the case e = Ly/L^ = 0.2. We are assuming, for this discussion that 
the two components of the gradients are equal 

dx dy 2 dz ^ ^ 

Otherwise each curve in fig. 6 will have to be multiplied by the appropriate gradient and 
the total effect will be altered. 

We see that around the zero crossing the contributions of the two dimensions contribute 
with opposite sign so the symmetry (equation 9) does not hold and if we wish to operate 
near the zero crossing the result will depend on and separately and not on ^ 
except in the case of cyhndrical symmetry ^ = ^ = \ ^ . 

Finally fig. 7 shows the comparison of theory and numerical simulation for T=300mK. 



IV. DISCUSSION 

The calculations of the velocity correlation function (vcf) presented here and in [4], start 
by following a single trajectory with the collisions only damping the amplitude of the vcf 
and not changing the velocity components, v^iVy as the particle moves along this trajectory. 
The simulations, on the other hand follow a particle as it is deflected to another (randomly 
chosen) trajectory by the collisions. Essentially the theory follows the particles that haven't 
collided while the simulations follow those that have. 

Each collision, while conserving the kinetic energy (magnitude of velocity), will result in 
a change of direction of the motion and hence of the x and y components of velocity so the 
correlation function resulting from the simulation of a single trajectory will be quite different 
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from that considered by the theoretical calculation. However, for each collision that takes 
a particle from trajectory (1) to trajectory (2) there should be a collision leading to the 
reverse transition according to detailed balance and when we average over all trajectories 
the results are seen to be the same. 

In addition while we have neglected the motion parallel to the E field we see that collisions 
will alter the velocity component in this direction, and, as a result of conservation of energy 
in the collisions, will thus alter the velocity in the perpendicular plane. However this effect 
will be cancelled when we average over all trajectories in the perpendicular plane as we have 
done here and in [4]. (We note that due to the heavy mass and slow He^ velocity, Baym 
and Ebncr [9] conclude that the phonon scattering on He^ is predominantly elastic.) 

We have presented the general solution for the frequency shift linear in E, for the case of 
a rectangular box with specular reflecting walls. The effects of non-specular wall reflections 
are expected to be small for the case of heavy damping of interest with respect to the co- 
magnetometers. For UCN the non-specular wall collisions are expected to be the major 
source of damping but the effect is expected to be small. This will be discussed in a 
subsequent work. 

V. FIGURE CAPTIONS 

Fig. 1) Normalized frequency shift vs normalized frequency, cu' — cUoLx/v, for the single 
dimension contribution of a single velocity. 

Fig. 2) One dimensional contribution to the normalized velocity averaged frequency shift 
vs. reduced frequency cj^: = uJoLx/ I3{T), for Temperatures T=0.1, 0.2 0.3 and 0.4K, using 
the temperature dependent mean free path for He^ in He^ 

Fig. 3) One dimensional contribution to normalized velocity averaged frequency shift 
with normalized frequency {uj* — ujLi/ (5{T)) as a parameter. 

Fig. 4) Normalized frequency shift for a rectangular box with e — Ly/Lr^ — 0.2 normalized 
to L^. vs frequency normalized to L^/ (3 (T) for various temperatures. 

Fig. 5) Frequency shift in a rectangular box with e = Ly/L^ — 0.2 normalized to vs 
temperature, K, with normalized frequency as a parameter. 

Fig. 6) Contribution of the long dimension, x, (red), short dimension, y, (blue) and 
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combned result (violet) for the frequency shift normalized to long dimension (x),vs, frequency 
normalized to the longdimension, col — u* — u;oLx/P{T) for a temperature of 0.4K and 
s = Ly/L^ — 0.2 The results of the theory are shown along with those derived by numerical 
simulations. 

Fig. 7) Contribution of the long dimension, x, (green), short dimension, y, (red) and 
combned result (violet) for the frequency shift normalized to long dimension (x),vs, frequency 
normalized to the longdimension, to* — uJoLx/l3{T) for a temperature of 0.3K and e — 
Ly/Lx — 0.2 The results of the theory are shown along with those derived by numerical 
simulations. 
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